Motion of a test particle in the transverse space of Dp-branes 
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We investigate the motion of a test particle in higher dimensions due to the presence of 
extended sources like Dp-branes by studying the motion in the transverse space of the brane. 
This is contrasted with the motion of a point particle in the Schwarzschild background in higher 
dimensions. Since Dp-branes are specific to 10-dimensional space time and exact solutions of 
geodesic equations for this particular space time has not been possible so far for the Schwarzschild 
background, we focus here to find the leading order solution of the geodesic equation (for motion 
of light rays). This enables us to compute the bending of light in both the backgrounds. We show 
that contrary to the well known result of no noncircular bound orbits for a massive particle, in 
Schwarzschild background, for d > 5, the Dp-brane background does allow bound elliptic motion 
only for p = 6 and the perihelion of the ellipse regresses instead of advancement. We also find that 
circular orbits for photon are allowed only for p < 3. 
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I. INTRODUCTION 

In classical theory of general relativity we study the gravitational effects of massive point particles (since far away 
from the source the gravitating object can be thought of as a point particle) and there are various tests of general 
relativity based on the analysis of the motion of point test particles (both massive and massless) in such a background. 
For a variety of reasons there has been substantial progress, over the years, in understanding general relativity in 
higher dimensional space-times. In fact, string theory is one of the leading candidates for the ultraviolet complete 
theory of gravity and the consistency of the super string theory requires that our universe should have nine spatial 
dimensions in contrast to the three spatial dimensions of the presently observed universe. This discrepancy is best 
resolved by the proposal of compactification which says that the extra six spatial directions are compact and the size 
of this compact space is so small that present experiments cannot possibly probe it directly. On the other hand, it is 
quite likely that the low energy theory would capture some of the effects of compactification from higher dimensions 
which can possibly be tested at LHC. From the analysis of a toy model of compactification from five to four dimensions, 
it is becoming overwhelmingly clear that in such a scenario, not only will the standard model particles be duplicated 
but the Kaluza-Klein modes of some of these particles will also have clear signals in experiments at LHC pj . While 
we await the confirmation (or absence) of such signals from the accelerator experiments, it is worth analyzing large 
scale tests of higher dimensional gravity. In fact, there already exist substantial literature on (higher dimensional) 
Kaluza-Klein gravity Q ■ This theory can be tested by analyzing the motion of test particles in the background of 
a static, spherically symmetric mass like the sun. The generalization of Birkhoff's theorem to higher dimensional 
compactified theories [3], 4] allows us to construct static solutions for the metric and in the five dimensional case 
most of the attention has been focussed thus far on the solitonic metric [H] , [f| , Q which satisfies the five dimensional 
vacuum field equations and which reduces to the standard Schwarzschild solution on a four dimensional hyper-surface 
without explicit dependence on the fifth dimension. 

The motion of test particles in the gravitational background of the soliton has been studied in connection with the 
well known classical tests of general relativity, e.g. gravitational redshift, light deflection, perihelion advance, time 
delay etc 0,0, [HJ- For the five dimensional case, the soliton metric is of the form 

ds 2 = A a dt 2 - A- a - b dr 2 - A 1 - Q - b r 2 (rf6» 2 + sin 2 Odtf) - A b dy 2 , (1) 

where y denotes the coordinate of the fifth dimension and A(r) = 1 — 2M jr. Here, M is the parameter (including the 
gravitational constant) related to the mass of the soliton taken to be at the origin of the coordinate system and the 
exponents a, b are constant parameters satisfying a 2 + ab + b 2 = 1. Thus the metric has two independent parameters. 
Note that we can recover the standard four dimensional Schwarzschild metric on the hyper-surface (y — constant) in 
the limit a = 1, b = 0. Hence b can be considered as the meaningful free parameter of the theory in the sense that a 
non-zero value of b signals a departure from the Schwarzschild geometry of a massive point gravitating source. 
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While the above analysis is motivated from the point of view of Kaluza-Klein gravity in five dimensions, we will 
be focussing in the following on the case of (super) string theory which is consistent in ten dimensions. This theory 
contains point like excitations in its perturbative spectrum as well as extended but stable objects called Dp-branes 
in its non-perturbative spectrum [Tlj. The Dp-branes define hypersurfaces in the ten dimensional space-time with 
p-spatial and one time like coordinates. Therefore, they are p + 1 dimensional objects with 9 — p transverse spatial 
directions and in this theory it will be of interest to study the gravitational effects, not only of massive point like 
sources but also of massive branes i.e. extended objects, on probe particles moving in the space transverse to the 
brane. This would allow us to compare and distinguish the effects of point like sources from those of sources which 
are extended in nature e.g. strings, membranes or other branes. 

Therefore, we are going to study the f O-dimensional curved space-time metric obtained in string theory describing 
such objects. Let us consider N coincident Dp-branes described by the metric (c = f ) [l2j 

ds 2 = dr 2 = ff-i (^dt 2 - X>z 1 ) 2 ^ - ffs (dr 2 + r 2 dQ 2 _ p ) , (2) 

where r denotes the radial coordinate and dQ\_ p the angular element of the (9 — p)-dimensional transverse space. 
The harmonic function H(r) is given by 

*(r) = 1 + nfe, (3) 

where Q p = d p Ng s £ 7 s ~~ p with g s ,£ s denoting respectively the string coupling and string length scale and d p — 
2 5 -P7r( 5 -p)/ 2 r((7 — p)/2). For all practical purposes we can assume that Q p ~ ( 7 S ~ P - We note here that 

H (r) 1, H(r)^A-^ : (4) 

where large and small r are defined with respect to the string length l s . The metric (|2|) satisfies the Einstein equation 
derived from the 10- dimensional gravitational action in the presence of a Maxwell term with a field strength of rank 
p + 2 (for notations and other details see [l2| ). We will study the behavior of a point test particle moving under the 
influence of this gravitational background in the space transverse to the p-brane. Let us note that the form of the 
line element given in ([2]) is clearly different in character from the Schwarzschild line element (to be discussed in next 
section) as well as the solitonic case in (JTJ, namely, the function H 2 multiplies both the radial and the angular line 
elements for the case of the branes which is to be contrasted with the other cases. Thus, the analysis we will carry 
out in this paper will have a different character from the five dimensional case described above. 

We point out here that there has already been an extensive study of the motion of probe branes in the transverse 
space of stacks of static branes of various kinds (see for example 03, [ll], |Hj], pj], UJ%, (Tljlll which have also 
been used in cosmological applications [2(|, [2l|, (22[, [23l] . [24j . Also, in mirage cosmology |25j.[2q. the motion of 
a D-brane is considered in a gravitational background. On the contrary, it is the analysis of the motion of probe 
particles which is the goal of the present investigation. In a sense, since higher dimensional gravitational theories do 
admit extended objects as gravitational sources, it is quite natural to investigate the motion of a probe point particle 
in the transverse space of branes (i.e. in the gravitational field produced by extended objects like branes) so as to 
contrast the results with those obtained from the motion of a point particle in the gravitational field produced by 
point sources. 

Before we proceed further, it is useful to discuss briefly about the distance scales at which the gravitational effects 
of these branes become important. We recall that in a space-time of arbitrary dimension d, the gravitational potential 
created by a point particle of mass M at a distance r is given by 

where denotes Newton's gravitational constant in (i-dimension. Working in natural units where V is dimensionless, 
we define a characteristic gravitational length scale L given by 

L d-3 = G (rf) M; ( 6 ) 

so that we can write V = — (L/r) . This shows that gravitational effects are weak for r » L but become important 
at scales of order L. Note that for a point gravitational source, L coincides with the Schwarzschild radius of a black 
hole of mass M in d-dimensional space-time (up to factors of order unity). 
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When we have a stack of N Dp-branes wrapped around a (spatial) p-dimensional compact volume V p , the effective 
mass of the branes is given by 

M = NTpVp ~ ^fi , (7) 

where T p is the tension of the Dp-brane and a' corresponds to the slope of the Regge trajectory (which is related to 
the string tension or the square of the string length scale £ s ). The dimensionally reduced space-time in this case is 
(d — p)-dimensional. In this space-time, the branes can be thought of as a point source of mass M. It follows from 
([51) that the characteristic size L of such a system (where p dimensions have been compactified) is given by [27j j 

L *-,-3 = G (*-P) M= g?M = G ^ N . (8 ) 

Using the relation between the Newton's constant in (^-dimensions and the string coupling 

G^~gUV^) d - 2 , (9) 

in © wc obtain 

L d-p-3 _ g a N(y/tf)*-P-S. (10) 

Note from (fTU)) that for d = 10 and for p < 6, we have L — > as g s N — > 0, i.e. the gravitational effects of the 
branes vanish. Thus, we need g s — » and a sufficiently large N such that g s N is finite for the validity of our 
analysis. We also note that when L < \foH , there does not exist any scale where the gravitational effects play any 
significant role. Besides, the fact that L in (ITU)) is independent of the compact volume V p , clearly shows that the 
concept of a characteristic length scale L is still relevant for the system of branes having infinite extensions in their 
tangential directions. In our analysis, we consider only the case where the probe particle is farther away from the 
branes compared to this characteristic length scale. 

The paper is organized as follows. To compare the results of the motion of a point particle in the transverse space 
of p-branes with that in the gravitational potential of a point mass in arbitrary d-dimensional space-time, we discuss 
the results for the latter case first in the next section. In section III, we study the motion of a point particle in the 
transverse space of p-branes. Last section is devoted to a brief conclusion and a discussion of the results. 

II. GENERAL RELATIVITY IN d-DIMENSIONS 

In this section we work out the motion of a point particle (both massive and massless) in the gravitational back- 
ground produced by a point like massive object. Let us note that there has been discussions, in the in the past, on 
solving geodesic equation for probe particles in Schwarzschild back ground. For example, in a recent paper [2g|, the 
geodesic equation for a spherically symmetric space time has been exactly solved. The solutions for d = 4, 5, 7 are 
presented in terms of elliptic functions and for d = 6,9, 11 are found in terms of hyperelliptic functions. However, 
for space time dimensions, 8, 10, > 12 , analytic solutions have not been possible. Since Dp-branes are specific to 
ten dimensional space time and we intend to compare various results of L>p-brane background with Schwarzschild 
background, we find in this section the leading order solutions (for motion of light rays) for arbitrary space time 
dimensions. Thus, let us consider a spherically symmetric gravitating object of mass M in c?-dimensions (where 
d > 4) located at r = 0. The gravitational effects of this object outside is described by the Schwarzschild metric in a 
rf-dimensional curved space-time [29l | given by 

ds 2 = dr 2 = h(r)dt 2 - J^jdr 2 - r 2 dfl 2 d _ 2 , (11) 

«0-(i-^). (12) 

Furthermore, d£l 2 d 2 in (1111) represents the spherically symmetric angular element and the interesting nontrivial 
components of the metric tensor follow to be 

90Q = h(r), gu = g„ = -— -r-. (13) 

h{r) 



where we have assumed c = 1 and 
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This metric describes the static and spherically symmetric vacuum solutions of the <i-dimensional Einstein's equations. 
It is easy to see from (fTIZj) that the Schwarzschild radius in this case is given by 

r { s d) = (2G^m)^~ 3 . (14) 

Let us note that if we compactify d — 4 spatial dimensions, as we have seen in ([8]). we can relate the d-dimensional 
Newton constant to the four dimensional one by the relation 

G (4 » = (15) 

Therefore, if the radius of compactification is very small, it follows from (fT5|) that for a given G^\ the higher di- 
mensional Newton's constant will be smaller, G^ <C G^ 4 \ Correspondingly, upon compactification the Schwarzschild 
radius (for the same mass) will also be smaller in higher dimensions, Tg <C ri 4 \ 



A. Motion of massive particles 



To study the motion of a massive point particle of unit mass in such a background, let us consider the Lagrangian 

L = g^x", (16) 

where the metric tensor is given by the Schwarzschild line element (|11|) and the "dots" denote derivatives with respect 
to the proper time (proper length) r. Since the metric in (fTTj) is static, x° — t is a cyclic variable and we expect its 
conjugate momentum to be conserved. In fact, the Euler-Lagrange equation for t gives 

d dL 

— r = 0, or, h(r)t — k — constant, (17) 

ar at 

which can be written equivalently as 

With this brief general introduction, let us study various aspects of a point particle motion in such a background. 
Let us make the following geometric argument to determine the number of coordinates that need to be considered for 
such a motion in general. We note that it takes three points to define a two-dimensional plane. For the first point, 
we can take the gravitating object assumed to be at r — 0. The second point will be our point test particle of unit 
mass which feels the gravitational force at coordinate r(t). We can choose the third point to be the same point mass 
only an infinitesimal time later at r(t + dt). These three points define a two-dimensional plane, and the gravitational 
force will never take us outside of this plane (because the gravitational force is central and pulls the outer point only 
towards the center radially). This can be understood from the point of view of angular momentum conservation as 
well. As a result, we can rotate our coordinates so that the metric and the Lagrangian have the simple two spatial 
dimensional form 

dr 2 = h(r) dt 2 - -^-dr 2 -r 2 d(f> 2 , L = h(r) i 2 - — ^f 2 - r 2 2 , (19) 
h(r) h(r) 

where h(r) is given in (pT2|) . We note from the Lagrangian L in (flT))) that not only is t a cyclic variable (in general 
for the Schwarzschild metric) as described in (|17p . but the angular variable <f> is also a cyclic variable. Therefore, the 
corresponding conjugate momenta are conserved and the i, 4> equations of motion lead to 

i = W) = ~ W*' (20) 



where k, £ are constants. We recognize £ to correspond to the angular momentum associated with the particle motion 
and at this point the two constants appear to be arbitrary. 
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1. Radial free fall 

Let us begin by studying the motion of a radially (vertically) free falling point particle of unit mass. In this case 
d<fi — (corresponding to the case I = 0), and dividing the line element in (fTTT) by dr 2 we obtain 

l = h(r)P-j^r 2 = J-^(fc 2 -f 2 ), or, r 2 = fc 2 - ft(r), (21) 

where we have used (p~8|) (or equivalently ([20|0 . If we assume that the particle falls initially from rest from a point 
ro, we have r = 0, we determine from (|2T|) 



P = Mr ), r 2 = fc(r„) - ft(r) - 2G«M ( -Lr - -L ) . ( 22 ) 



Taking the r derivative of (l2"2")l we obtain 



2rr = -(d-3) r, 
(d-3)GWM 



(23) 



which is the analogue of Newton's equation in d-dimensions. Clearly, in higher dimensions the gravitational force falls 
off much faster for large r and the attraction is much stronger for smaller r (compared to 4 dimensions) . 

From (|2"21 we note that the proper time that a particle will take to come to a radial coordinate r (at r) is given by 

t r tq _l 

T= [ dT >=[ d C=[ d r' ( - V • (24) 

J J r' J \2G^M(4- 3 - (r') rf - 3 )/ 

ro r 

Here the sign of r (equivalently the limits of the integration) are chosen such that r increases as the radial coordinate 
r decreases. Since the integrand is well behaved, the final integration limit can actually be extended all the way to 
the Schwarzschild radius (assuming that it lies outside the gravitating mass) and the proper time taken to reach there 
is easily seen to be finite. 

Let us also note that the coordinate speed of the particle is given by 

V( ' ~ dt~ i~ / x 2G(")M ^ (r r) d - 3 j ' [ ^ 

which shows that the coordinate speed of the particle vanishes at the Schwarzschild radius (as well as at r = ro 
initially) . Since the speed is increasing initially, it must reach a maximum speed at some point before decreasing and 
the radial coordinate where the coordinate speed is maximum is obtained from 

dv(r) n ( 3G^M4~ 3 
— — = 0, or, rw = ^ 



where the speed has the value 



2 / 2G^M\ 2 l2 



All these results can be easily seen to reduce to the standard results when d = 4. From (|26p and (|27p we see that for 
ro > r s (c = 1) 



(qG^m) ^ = (3)^ r s , w = V (r max ) = * (28) 



3^3' 

This shows that, in this case, the maximum speed is reached very close to the Schwarzschild radius. In fact, 1.7rg < 
r m ax < as 10 < d < 4 and for ro 3> r s the maximum speed achieved by the test particle is independent of the 

number of space-time dimensions. 
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2. Circular orbit 



We will now discuss more general solutions of the particle motion. First, let us consider circular motion of the 
particle. To begin with let us note that classical Newtonian gravity does not permit bound circular orbits in higher 
dimensions d > 4. To see this, let us recall that the Lagrangian for the Newtonian motion in d dimensions is given by 
(the test particle has unit mass) 

-i((£)"+4(S)>f£ 

and the energy associated with the particle motion is obtained to be 

1 / fdr N \ 2 , fd(j>\ 2 \ G^M , s 

E =Hw +r "{i) )- — <0 - (30) 

for bound motion. The ^-equation leads to the conservation law 

Tt = 4". (3D 

dt r 2 N 

and the radial equation for circular motion (^f- = = < ^,!j v ) leads to 

'd<j>\ 2 (d-3)G^M 



dt 



r d-2 

I TV 



(32) 



Equation (|32[) simply expresses the fact that, for circular motion, the centrifugal force must balance the gravitational 
attraction. Using pip we can determine the radius of a possible circular orbit r N from (|32p to satisfy 



. <£HpL. (33, 

Clearly, there are two distinct cases to consider. For d = 5, (|3"3")l leads to the constraint 

£ 2 = 2G (5) M, (34) 

for which we determine from (1301) that 



£7 = 0. (35) 



For d > 5, (l33l) determines 



(36) 



which leads to the energy associated with motion to be 



£ 2 d-5 

E=^r— 3 >0, d>5. (37) 

In either case we note that this does not correspond to bound motion for which E < 0. 

This can be understood more physically as follows. Let us define an effective potential associated with the particle 
motion as (see (|5D])) 

V«{r») = - -PS-, (38) 

where the first term represents the centrifugal potential which is independent of the space-time dimensions while the 
second term representing gravitational attraction does depend on d. When d — 4 the positive centrifugal potential 
dominates the gravitational attraction for small values of r N while for large r N it is the other way around. This leads 
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to a stable minimum of the potential at some negative value (for (|30l) to be true with r N constant) for a finite r N 
where the centrifugal force can equal the gravitational attraction. For d = 5, both the terms in the effective potential 
have the same -3- behavior and if £ 2 — 2G^M, the effective potential vanishes and is incompatible with (|30l) . For 
d > 5, the gravitational attraction dominates the positive centrifugal attraction for small values of r N while it is 
the other way around for large r N . As a result, the extremum of the effective potential in this case is an unstable 
maximum and occurs at a positive value of the potential (centrifugal potential is positive) which is inconsistent with 
(f30|) . Alternatively we note that the extremum of the effective potential is derived from 

, . . 1 / , (d-3)G^M\ , . 

V: s (r N ) = -~jU 2 - = 0, (39) 

' TV \ I AT / 



d=5: £ 2 -2G^M = 0, d > 5 : r N = ( {d - ^f^ ) ^ . (40) 



which determines 



This can be compared with the earlier results in (|34[) and ()36[) . Furthermore, the second derivative of the potential 
at the extremum gives 

d=5: V^rv) = 0, d > 5 : V&(r N ) = ^ (5 - d) < 0, (41) 

T 

implying that the extremum is an unstable maximum (for d > 5) which cannot support bound state motion. 

On the other hand, inclusion of relativistic effects does allow for stable circular motion. For example, the r-equation 
following from (Til)]) leads to (r = r = 0) 

rr-<±JjgM.e, (42) 

and this equation (compare with (|32j0 clearly shows that the two constants defined in (|20|) are related. Explicitly we 
have 

f = (d-3)GWM fe 2 m 

r d - 5 (. 2G wm \ 2 ' 

I 1 r d-3 j 

Using (|42p in the relation following from the line element for circular motion, namely, 

1 = h(r)i 2 -r 2 ^ 2 , (44) 
as well as the relations shows that (l4"2j) leads to 

(d-l)G^M\ k 2 



_> 



(45) 



This relation is consistent only for 



2 

showing the possibility of circular motion. In fact, (|45[) determines the value of the constant (for a fixed r) to be 

/ -i 2G W M \ 

I 1 ,.£(-3 I 

k = ; 

, _ (d-l)G< d )M 

y *■ r d~3 

To understand circular motion better, let us use P2^1 to determine 



r>(M^ r « (46 ) 



(47) 



'(d-3)G( d )M 

dt ~ 1 



<ft= V (d-3)GWM I**' (48) 
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where we have assumed r is a constant for circular motion. Equation (|48|) determines the (coordinate) period of the 
orbit of such a motion to be 



T = At = 2tt —7- — . (49) 

This can be thought of as the generalization of Kepler's law to d-dimensions, namely, the square of the period of the 
orbit is proportional to the radius of the orbit raised to (d — 1) power which corresponds to the number of (transverse) 
spatial dimensions in d space-time dimensions and this reduces to the standard result when d = 4. (One should keep 
in mind though that these are only coordinate parameters.) 

To determine the proper period associated with such a motion, we note from the first relation in (I20[) that we can 
write 

/ 2G&M\ At 
^={l-^^) T , (50) 

and using the value of k from (|4"7|) leads the proper period ([50)1 to correspond to 

^(i-^fca ) A '- (51) 

We note that the proper period vanishes for 

'V-^r^, (52) 

which would correspond to the circular orbit for the photon in d-dimensions (which has a vanishing proper time) and 
this coincides with the well known result when d = 4. 



3. General motion 

If we do not assume the particle to be moving in a circular orbit, the radius will not be a constant so that r =/= 0, r =/= 0. 
In this case, the r equation is different from (j42|) and, consequently, the two constants k, £ defined by (|20|) are not 
directly related. In fact, the r equation can be easily obtained from the line element in (fT9]) which leads to 



(i- 2 ^) / i 2GW>M \ 2 P r2 f 2GWM 

r d-3 J ^ 2 r \ r d-3 

p \f 1) + -^r) [1 —) , (53) 

where we have used the two relations in ([20|) and £ is assumed to be nonzero. [I = corresponds to the case of radial 
fall which we have discussed earlier. For I and k defined by (|4"3")) the right hand side of (1531 vanishes leading to a 
constant r, namely, a circular motion.) Defining u = -, this equation leads to 

fdu\ 2 2G^Mu d ' 3 „/ r] ,\ (fc 2 -l) , x 

Upon differentiation, this gives two possible equations 

du „ d 2 u (d-3)G^M d 4 . , ,,_ f(ft , # d 2 

— =0, or, — s = i ^ ti"-ti + ((f-l)G (£ ' ) Mu". (55) 



The first equation in (|55p corresponds to the circular motion which we have already discussed. The second equation, 
on the other hand, describes the more general motion of a massive particle. In d — 4, we know that the first two 
terms on the right hand side lead to elliptic orbits whereas the last term corresponds to relativistic correction to these 
orbits resulting in the perihelion advance. Correspondingly, the equation 

(d-3)G^M 

#2 - u N + (2 u N , (56) 
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can be thought of as representing the classical Kepler's equation in o?-dimensions for a massive particle moving in the 
gravitational potential ([5]) . The last term in (|55[) would then represent the relativistic correction to this equation. 
The Newtonian relation for bound motion analogous to (|54|) (or (|30|) ) has the form 

Equation (155)) is solved in closed form for d = 4, 5, 6, 7, 9, 11 and could not be solved for other values of d psj . However, 
the analysis of the stability of any bound orbits can be done easily from (|57p (see (|3T?]) - PTj) ). Since d = 4 is well 
studied, let us restrict to d > 5. In (1571) we can identify 

Kff(M«) = Ujv ^2 ' ( 58 ) 

which coincides with (1381) (up to a multiplicative constant) and, therefore, the stability analysis of the motion goes 
exactly as in (|39[) - (|41[) showing that no bound motion is possible. Intuitively this can be understood from the fact 
that for r»r s (meaningful for physical motion) the attractive gravitational potential falls off much faster for d > 5 
than the repulsive centrifugal potential leading to unbounded motion. 

On the other hand, including the relativistic correction, we can define the effective potential from ([M]) to correspond 

to 

V cS («) = - 2G(rf ^' 3 + v 2 (l - 2G^Mu d - 3 ) . (59) 
The extremum of this potential is obtained from 

V: B (u) =2u(l- (rf ~ 3 ]f rf)M u d - 5 - (d - l)G^Mu d -^j = 0. (60) 

This equation has one vanishing root and only one positive root uq (which is physical) for d > 5 (this follows easily 
from Descartes' rule of signs). The second derivative of the potential at the positive root leads to 

V&(uo) = -2 ((d - 5) + 2(d - l)G (d) Mu^ 3 ) < 0, d > 5, (61) 



leading to instability of motion. Once again, we see that the relativistic correction in (|59j) only adds an attractive 
potential which falls off even faster than the classical gravitational potential and thereby does not help in stabilizing 
the motion. Thus, unlike the circular motion, including relativistic correction does not stabilize the general motion 
and this was the main thrust of the analysis in [29| . 



B. Motion of light rays 

Finally, we discuss the motion of photons in the higher dimensional gravitational field. As we have already seen in 
(|52[) . photons can move in circular orbits of radius (as we have emphasized earlier, our analysis is for d > 4) 

r = ^ rW> = ((d - 1)G^M) ^ . (62) 

Normally, such a radius lies inside a star and, therefore, is not physical. There is no other bound motion possible 
for the photons. Therefore, in this section we will discuss the phenomenon of bending of light in the vicinity of a 
gravitating star in d dimensions. 

Since photons are massless particles, the proper time associated with them vanishes and we have dr 2 = 0. Conse- 
quently, we cannot label the photon trajectory with r and have to introduce an affine parameter A for this purpose. 
Equations (|20[) continue to hold 
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where the dots now refer to derivatives with respect to A. Only the radial equation modifies which can be obtained 
easily from the relation for the line element. Since dr 2 = it follows that 

dr\ 2 ( 2GWM\ , fcV . . 

In terms of the variable u — -, this can be written as 

^V + (l-2GWMu rf -V = ^, (65) 

which can be compared with ([54]) . Taking derivative of this equation with respect to u we obtain 

du ( d 2 - 



d<p 

which leads to two possibilities 



= 0, (66) 



S= ' ^+u-(d-l)G^Mu d - 2 =0. (67) 
dip dtp z 

The first of these describes the circular orbit that we have already discussed. Therefore, let us analyze the second 
equation 

^ + u = (d - l)G^Mu d - 2 = eu d - 2 . (68) 
d<p z 

where we have identified the relativistic correction as 

(d - l)G {d) Mu d - 2 = eu d - 2 , (69) 

and recognize that it represents a very small correction compared to the centrifugal force for r 3> ri ■ 



The solution to equation (|68|) without the relativistic correction is that of a straight line which can be written as 

u(4>) = u min sin 0, (70) 

where we have chosen the closest distance of approach of the photon to the star r m ; n = to occur at <j> = S . Using 
this we can write down the leading solution (in powers of e) to (f68|) of the form 

u(4>) = M min sin0 + eu\(4>). (71) 

Substituting this back into (|68l) we obtain to linear order in e 



d ui 



ui = t4 n 3 frin''-V (72) 



The solution of this equation in even dimensions can be written in the form 

d 

2 

Ml (0) = £a„sin^ 2 "0, (73) 

n=l 

with the constant coefficients given by 

w d ~ 2 (d — 2n) d 

ai= - (d-l)(5-3) ' a " +1= (d-2n-3) a - n = l,2,...,--l. (74) 

Therefore, the complete solution to the leading order in e can be written in even dimensions as 

d 
2 

u(<p) = u min sin<?(> + (d - l)G {d) M "n sin rf - 2 ™ 0. (75) 

71=1 
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This is easily seen to reduce to the well studied results when d = 4. 

Let us next note that in the absence the relativistic correction, the solution ([70)1 shows that the light ray would 
come in from infinity at an angle = and go out to infinity at <p = n. In the presence of the relativistic correction 
we expect the light ray to come in from infinity at an angle <j) = —5 and go out to infinity at (j) = tt + 8 (by symmetry) 
where we expect the angle 8 to be small. Using this in (|75|) we obtain 



= u min sm{-5) + (d-l)G {d} Mj2a n sm d - 2n (-5). (76) 

n=l 

Recognizing that 8 is a small angle, we determine 



\ 

n=l 

The bending of the light path, therefore, is given by 



We note here that the overall negative sign reflects the fact that the product on the right hand side is negative. 
In odd dimensions (recall that d > 5), the solution of (1721) can be obtained in a similar fashion and has the form 

Ul ((f>) =/3<f> cos 4>+ ^ a n sm d - 2n (t), (79) 

71=1 

where [|] denotes the maximum integer part of i and the coefficients satisfy 

u d ~ 2 (d — 2n) d 

ai = - (d-l)(d-3Y a ^- (d-2n-3) a - ^ = 3Q iiM' n = 1,2,..-, [-]-2. (80) 

In this case, we can also carry out the calculation for the photon deflection in an analogous manner. The equations 
force that the light ray comes in from infinity at <j> = and goes out at (j) = tt + 8 and the leading value of the total 
bending is determined to be 



III. MOTION IN TRANSVERSE SPACE OF Dp-BRANES 



In the previous section we systematically analyzed the motion of a point particle in a Schwarzschild background 
in higher dimensions and discussed various consequences associated with it. In this section we will carry out the 
corresponding analysis for the motion of a point particle in the transverse space of a stack of Dp-branes in ten 
dimensions (d = 10). The gravitational field produced by the stack of branes is given in ^ and, as before, we can 
show that the motion of the test particle will be planar so that the invariant length of the probe particle takes the 
form (c = 1) 

ds 2 = dr 2 = H-^{r)dt 2 - H^{r) (dr 2 + r 2 d(j) 2 ) , (82) 

where the harmonic function H{r) is defined in (|3]). The nontrivial components of the metric tensor are obtained 
from (|82|) to correspond to 



.9oo = H 2 (r), g n = g rr = 



-H*{r), g 22 =g 4>4> ^-r 2 H2. 



(83) 



12 



A. Motion of massive particles 

Let us assume that the test particle has a unit mass so that the dynamics of the particle is described by the 
Lagrangian (TTB")) 

L = g^± v , (84) 

where i M = and the components of the metric tensor g^ v are given in (|83|) . As before we note that t and <f> are 
cyclic variables so that the corresponding conjugate momenta are conserved. From the t, (f> equations, it follows that 
(compare with (j2T)]) ) 

i = — = kHi(r), 4> = — — = _.#-* ( r ), (85) 
5oo 94><i> T 

where k and £ are constants. With this let us now analyze the different types of motion the (massive) point particle 
can have in such a gravitational background. 

1. Radial free fall 

Let us start with the radial free fall of the point particle in this gravitational background. In this case we have 
I = since only the radial coordinate changes (d<f> — 0) and from (|52"]l we obtain (using (|85p ) 

1 = m(r){k 2 - r 2 ), =*> f 2 = k 2 - H~?{r). (86) 
If we assume that the particle falls from rest from a point ro, we have r\ ro =0. As a result, we see from (|86[) that 

k 2 = H-i(r ), r 2 =H-i(ro)-H-i(r). (87) 
Taking the r derivative of the second equation in (|8T[) we obtain 



which can be thought of as the analogue of Newton's equation in such a background. This equation can be compared 
with ([2^)1 and it is clear that as p increases, the gravitational force due to the branes falls off more slowly for large 
r than that due to the Schwarzschild background (for d = 10). This is a reflection of the fact that the number of 
transverse spatial dimensions decreases for higher dimensional extended objects. For short distances (see (j4])), the 
gravitational force, in fact, decreases for p < 5, becomes a constant for p = 5 and grows slowly for p = 6. This should 
be contrasted with (|23p where the gravitational force increases quite strongly as r becomes small. 
The coordinate velocity can be obtained from (1571) and is given by 

v(r) = ^ = ~ i = -Hi(r )H-i(r) (fl-*(r„) - H^{r)) ' = -fT* (r )G(r) (1 - G(r))* , (89) 
where we have defined 

i 

G(r) = HHr )H'i(r) = (j^y) * ■ (90) 



- "(2 ~ 3G(r)) ^ ^) 3 « = 0, (91) 



It now follows that 

_ dr v " "" v \4:H(r )(l - G(r)) J dr 

for G(r max ) = |. Since the radial coordinate is decreasing with time, H(r) > H(ro) and such a solution is possible. 
The coordinate where the velocity will be a maximum is given by 



4Q P 4r^ p / V 5 
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and the magnitude of the (maximum) velocity at that point is given by 

Kw)| = ff-*(r )G(w)(l - G(r max ))^ = -jLfH(ro) = ^= fc 2 , (93) 

which can be compared with ([27)1 . However, note that for a given ro, the constant k 2 is different for the Schwarzschild 
and the brane backgrounds. Nonetheless, it is worth emphasizing that in either case, for large enough ro, the maximum 
speed a test particle can achieve in a radial fall is given by (see also and note that c = 1) 

= -£j=. (94) 

The other significant point to note from (|89[) is that unlike the Schwarzschild background where the point particle 
has vanishing coordinate velocity at the Schwarzschild radius (in addition to the initial point, see t|25|) and discussion 
there), here the coordinate velocity does not seem to vanish (besides the initial point) except at the singular point 
r = 0. 



2. Circular orbit 



Let us next consider a slightly more general motion of the massive test particle (of unit mass) , namely, the motion 
in a circular orbit around the gravitational source. The line element (I82[) . in general, leads to 



I = H-2(r)t 2 - Hz(r) (r z +r z <j) z ) , (95) 

where dots denote derivatives with respect to the proper time r. If we are interested in a circular orbit, then r = = r 
and (l9"5l) leads to 



1 = H-i(r)i 2 -r 2 Hi{r)4> 2 . (96) 
On the other hand, the radial equation from (fT6|) (or ([84)) ) leads to 

or or 

Equations (|96| and (|97l) show that the two constants k, £ are not independent in this case. Explicitly, using (|85l) in 
(|M| we obtain 



t = 5 (r) ( ri? 5 (r) - lj , (98) 
which can be compared with (|43[) . Equation (I97[) allows us to determine 



AH (r) + rn'yr) 

where a prime denotes a derivative with respect to r. Substituting this into (I96p and using the definition of k from 
(1551 we determine 



2H~i(r) V 2H(r)+rH'(r 
We note from equation (|99| that we can write 



l) J =< 4r '^- 3W >>S- <ioi) 

It is clear from (jlOlf) that the right hand side is positive for any value of r for p > 3 (recall that we are considering 
< p < 6, see discussion following (Q|) and, therefore, circular orbits are possible in general in these cases. When 
p < 3, the right hand side is positive only for 

r> (j!=d&y\ (102 , 
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and circular orbits are allowed only for these values of the radial coordinate. The period of the orbit can now be 
determined from (|10ip to be 

T = At = 2n { (7-p)Q p rt-P ( 4r? ~ P + {P ~ 3)C ^ ^ + Q ")) 2 ' (1 ° 3) 

which represents the generalization of Kepler's law for such a gravitational background. It does not have a particularly 
simple form. However, for large values of r we note that 

T^r 2 ?, (104) 

which can be compared with (J49j) (with d = 10) and shows that the square of the time period has a smaller power 
growth with radius as the brane dimension increases (compared with the Schwarzschild background) . In general, for 
any background and for large radial coordinates it seems that Kepler's law generalizes as the square of the period 
is proportional to r a where a denotes the number of transverse spatial dimensions associated with the gravitational 
source. Finally, we note that the proper period associated with the circular orbit can be determined to have the form 

It is interesting to note that the prefactor multiplying At in (|105|) vanishes for p < 3 precisely for (see also ()102j) ) 

(3-p)Q P 



(106) 

and would correspond to the radius for the circular orbit of a photon. However, for p > 3, the prefactor does not 
vanish and correspondingly, the photon cannot have a circular orbit for such Dp-branes. 



3. Perihelion advance 



Let us next analyze the general (noncircular orbial) bound motion of a massive point particle (of unit mass) in the 
background of a stack of Dp branes. From the equation for the line element (|95j) we have 

(£)' 

Using the relations in (|55|) . we can write (|107l) as 

Gl) 2 + r2 i Hh (r) ( fc2jff * (r) _ x ) = 0) (108) 

and this shows that for k,£ satisfying (|98|) . this leads to a circular motion which we are not considering. As before, 
let us introduce the variable u = - so that (|108[) can be written as 

(|) 2 + " 2 - r^~ v Q * u7 ~ p y = i- (io9) 

This can be the starting point of our analysis for the stability of motion as we have done in subsection III A 31 However, 
because of the square root, a general analysis of this relation becomes quite complicated. On the other hand, we note 
that for physical motion we expect Q p u 7 ~ p -C 1 and, in this case, we can approximate 

(1 + Q p u 7 -") * = 1 + - 9i u ^p + ^u 21 - 3 " + 0(Q*), (110) 

and study (|109p perturbatively. As we will see, any motion of the perihelion can be seen only at the third order in 
perturbation. 
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We note that if we substitute the leading order term (up to linear order in Q p ) from (| 1 10|) into (|109l) we obtain 

du\ 2 , _ (2fc 2 -l)Q 7 _ v (k 2 -l) 

u 2 t 2 • 1 1 

From (|llip we can define the effective potential in which the particle moves as 

V cS (u)=u"- {2k2 ~f Qp u^, (112) 
so that the extremum of the potential satisfies 

rrl I \ „ A (7- P )(2fc 2 -l)Q p 5 \ 

KsW = 2u I 1 — 2 u° p J = 0, (113) 

where a prime denotes a derivative with respect to u. We note from (|113l) that for k 2 > ~, Descartes rule of signs 
determines that the equation has only one positive root (physical solution) which is given by (we are considering 
P < 6) 

M ° = l(7- P )(2 fc2 -l)0j ' (1M) 

On the other hand, for k 2 < i, there is no positive root leading to the fact that bound motion is not possible in this 
case. Therefore, we restrict only to the case k 2 > i where a physical root of (|113[) is possible. To determine the 
stability of this solution, we note that 

v»( \ (7- P )(5-p)(2k 2 -l)Q p 5 „ p 

Kff( w o) = ^ --2(5-p). (115) 

This relation is very interesting because it shows that the stability of the motion depends on the type of branes that 
provide the gravitational background. For example, for p < 5, we see that V^(uq) < so that the motion is not stable 
and we cannot have any bound orbit. For p = 5, the second derivative vanishes so that motion is again unstable. Only 
for the case of p = 6, we have V^(uq) — 2 > so that stable motion is possible. In fact, this seems distinctly different 
from the case of a Schwarzschild background where, as we have seen in subsection III A 3[ no bound motion is possible 
for d > 5. However, in light of the discussions in ([56 |) -([6T |) . we note from (|112p that with increasing dimension of the 
extended gravitational source, the fall off of the gravitational potential becomes slower and overtakes the centrifugal 
potential for p = 6. This is the reason behind the stability of motion in the background of Dp-branes. 

Therefore, let us consider the particle motion for k 2 > ^ and p = 6. In this case, (jllip as well as the equation 
following from this (we ignore the circular motion) can be written as 

du\ 2 2 (2k 2 -l)Q 6 (k 2 -l) 

— I + U — U = J 2 > 

(116) 



) 2t 2 I 2 

d 2 u (2k 2 - 1)Q, 



,2 



4£ 2 



This has the same structure as the motion of a point particle in a Schwarzschild background in d — 4 and the solution 
satisfying the two relations in (|116|) is given by 



(2k 2 - 1)Q 



u(4>) = M 2^ 6 (1 + e x cos < 



16(fc 2 - \)l 2 
(2k 2 - 1) 2 0I 



ej = l+ ,o,9 1^2 . ( 117 ) 



and we note that the bound motion is given by an ellipse with eccentricity e\ < 1 for k < 1 (in fact, we should have, 
in addition, ^k^-iye}'^ < ^ ^ or P n y s i ca l motion). Therefore, we will restrict ourselves top = 6,l>fc 2 >i (including 
the more stringent restriction given above) and in such a case, the perihelion and the aphelion are given by 

rAH - (2k 2 - l)0e 1 ~ ei ' PH " (2fc 2 -l)0 6 l + er 1 * J 
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If we add the next order (second order in Qq) contribution from (jllpp . the two equations in (|116[) generalize to 

du\\ A Q§\ 2 (2k 2 - 1)Q 6 (fe 2 -l) 

)2 



U 2 J 2i 2 t 2 



u ( Q 2 \ (2fc 2 -l)Q ( 



which have the same structures as (|116l) except for a small correction in one of the coefficients. The solution satisfying 
both these equations is given by 




- M 2 

(-0) 



(2k 2 - \) 2 Q 2 \ 8£ 2 
and we see that the motion continues to be that of an ellipse with 

— (2fc 2 — 1)Q6 \ 8£ 2 j 1 - e ' V2 ~ rpH _ (2k 2 - l)Q e ^ 8i 2 Jl + e' 
On the other hand, if we include the next order (cubic in Qg) correction from (|110p . then (| 109[) takes the form 

d„V _(*»-!) , (2*>-l)Q a « f 1 _Q£) u 2_ eu 3 e= M- m (1 22) 



t 2 2l 2 \ 8£ 2 J ' 16^ 2 

This has exactly the same form (with some differences that we will elaborate on in the following) as that for a point 
particle in a Schwarzschild background in 4 dimensions. We note that the locations of the perihelion and the aphelion 
are obtained by setting = so that the right hand side leads to a cubic equation. If we denote the three roots as 
Mi, U2, i*3, then it follows from the form of the cubic equation that 

Ul +u 2 +u 3 = --(l-^], (123) 



c 



8£ 2 



which is a large negative constant. Since we know that Ui,U2 are finite positive quantities (see (I12ip ). it follows that 
U3 must be large and negative and the motion is bounded by u\ < u < Following the standard method, we can 
now determine from (|122[) 



dcf> 1 — f(u + u\ + tt 2 ) 
du ~ yj(u-u{)(u2 -u\ 

which can be integrated to give 

1 - |(u + Ul + u 2 ) 



(124) 



It follows, therefore, that 



|^1= / ^^= 2V " ' =v(l--(u l + U 2)). (125) 



ire . . irffl/1 ,1 



A4> = 2M -2n = - — («! + u 2 ) = - + - . (126) 



2 v 1 ' z/ 32^ 



n r 2 



Namely, at the end of a period, the perihelion does not come back to its old location, rather, there is regression 
in its position. This is interesting because in the standard Schwarzschild case, the perihelion advances at the end 
of a period. This difference can be understood from the fact that while in the Schwarzschild case, the relativistic 
gravitational interaction is attractive, here the cubic term in (|122p leads to a repulsive interaction which slows the 
particle and hence leads to a regression. Of course, the motion of the perihelion is extremely tiny as in the case of the 
Schwarzschild background in 4 dimensions. However, this qualitative difference in the behavior of a brane background 
from that of the Schwarzschild background is most certainly interesting and any possible test for the presence or 
absence of this would be quite remarkable. 
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B. Motion of light rays 

Let us next discuss the motion of light rays (photons) in the background of a stack of Dp branes. As we have 
discussed earlier, in this case dr 2 = and the trajectory has to be parametrized by an affine parameter, say A. As a 
result, the line element leads to 

tf-i(r)t--#l(r) i^+r 2 ) = 0, (127) 



where dots denote derivatives with respect to the affine parameter A. Using the definitions in f|85[) we can write this 
relation as 

dv \ k 2 

' + r2 - ^ R ( r ) = °- ( 128 ) 



Using the new variable r = — , this equation can be written as 



-^ + u 2 -^u^ = ^, (129) 



which can be compared with (|55|) . If we differentiate this relation with respect to cf>, we obtain 

y \„:^ p u e - p ) =o, (i3o) 



which leads to 



du ( d 2 u k 2 (7-p)Q p 



— -0 or — + u - J ^L—P}Rp^-p - (ml 
dS-" 1 OT ' d6 2+U 2i 2 U ~ U ' {6) 



The first equation in (I131[) describes any possible circular orbit for photons which we have already discussed in (|106l) . 
Therefore, let us look at the more general second equation in (I13ip . 
We note that the second equation in (|131[) can be written as 



,2 



+ u = eu 6 - p , (132) 



where we have identified e = - , which is a very small parameter. Comparing with ([68]) we see that (|132[) has 

the same structure with the identification d = 8 — p and the analysis in (|69 ]) -(|8T |) can be carried through completely 
in a similar manner. The deflection in the path of a light ray due to the presence of a stack of Dp branes can now be 
determined to be 

%S II (l + l=jhs), Peven, 

T-m ■ ( 133 ) 

3J S^ rf (1+53^), POdd 

n— 1 v / 



IV. CONCLUSION 



In this paper, we have systematically analyzed the motion of a point particle (massive and massless) in the grav- 
itational background of a stack of Dp-branes in 10 dimensions. Such extended objects exist in higher dimensional 
theories of gravity such as the string theory and, therefore, it is quite natural to study the effect of such gravitational 
sources on test particles and to compare and contrast them with the effects produced by a Schwarzschild background. 
For this purpose, we have also systematically analyzed the effects of the Schwarzschild background on test particles 
in higher dimensions (some of which may be available in a scattered form) and then carried out the similar analysis 
for the Dp-brane backgrounds. In both cases, we have studied the radial fall of a massive particle, the motion of a 
particle (massive and massless) in a circular orbit as well as the general motion of a massive particle. In addition, wc 
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have also studied the question of the deflection in the path of photons in such backgrounds and we summarize our 
results in the following. 

The analysis of the radial fall of a massive particle in the Schwarzschild background leads to the well known fact 
that the (coordinate) speed of the particle falling initially from rest at ro keeps on increasing up to a point and then 
starts to decrease until it vanishes at the Schwarzschild radius (in d dimensions). In the background of the Dp-branes, 
the particle falls from rest and gains speed up to a certain point and then the speed decreases, but vanishes only 
at the origin. It is interesting that independent of the number of dimensions of space-time and the nature of the 
gravitational source (namely, Schwarzschild or Dp-branes), the maximum speed that the particle can have (for large 
ro) is given by v max = 

In a higher dimensional Schwarzschild background, circular orbits for a massive point particle are possible if the 
radius of the orbit satisfies 

r>(<iZ») A r!f>, (134) 



where d denotes the number of space-time dimensions. For a background of Dp-branes, on the other hand, circular 
orbits are possible without any restriction on r for p > 3 while for p < 3, a circular orbit requires the radius to satisfy 



In all the cases, we have shown that for orbits with large r, a generalization of Kepler's law holds, namely, the square 
of the period is proportional to r Q where a denotes the number of transverse spatial dimensions associated with the 
given gravitational source. 

For a massive particle in a Schwarzschild background there is no noncircular bound orbits for d > 5 (so that the 
motion of the perihelion of an elliptic orbit is not a meaningful question) . This was the main result in the work of [29j 
and this can be understood qualitatively from the fact that for r » r Sl the attractive gravitational potential falls off 
much faster than the repulsive centrifugal potential leading to instability in motion. In contrast, the background of 
Dp-branes does allow bound elliptic motion for p = 6. The reason for this lies in the fact that as p increases, the fall 
off in the gravitational potential becomes slower and overtakes the repulsive centrifugal potential precisely for p = 6. 
Consequently, one can study questions such as the perihelion motion in this case. We find that when the relativistic 
corrections are included, the perihelion of the ellipse regresses (as opposed to the conventional advancement) and 
this can be traced to the fact that the gravitational potential resulting from the Dp-brane metric leads to a repulsive 
potential at this order contrasted with the attractive potential resulting from a Schwarzschild background. This is a 
significant difference between the two backgrounds and if this can be tested that would indeed be quite remarkable. 
We note that like the conventional perihelion motion, this effect is extremely tiny. 

For massless particles (photons), the Schwarzschild background permits a circular orbit for 

Mz^T^rW (136) 



which lies very close to the Schwarzschild radius in any space-time dimension. For the background of Dp-branes, on 
the other hand, a circular orbit for photons is allowed only for p < 3 at the radius 

(3-p)Q„ 



r=[^^) , (137) 

whereas there are no circular orbits for photons for branes with p > 3. Bending of light in the presence of these two 
backgrounds can also be calculated and are similar in nature. 
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